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We propose a scheme to perform a fundamental two-qubit 
gate between two trapped ions using ideas from atom inter- 
ferometry. As opposed to the scheme considered by J. I. Cirac 
and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995), it does not 
require laser cooling to the motional ground state. 



Quantum computation allows the development of 
polynomial-time algorithms for computational problems 
such a prime factoring, which have previously been 
viewed as intractable [y . This has motivated studies into 
the feasibility of actual implementation of quantum com- 
puters in physical systems 0. The task of designing a 
QC is equivalent to finding a physical realization of quan- 
tum gates between a set of qubits, where a qubit refers to 
a two-level system {|0), |1)}. Any operation can be de- 
composed into rotations on a single qubit and a universal 
two-bit gate, e.g. a gate defined by 

Ci2:|ei)|e2)^|ei)|(l + ei)©e2) (ei,2 = 0,l). (1) 

Implementation of a quantum computer requires precise 
control of the Hamiltonian operations and a high degree 
of coherence. Achieving the conditions for quantum com- 
putation is extremely demanding, and only a few systems 
have been identified as possible candidates to build small 
scale models in the lab Q-|4| . One of the most promising 
examples is a string of cold ions stored in a linear trap [^ . 
In this ion trap quantum computer qubits are stored in 
long-lived internal atomic ground states. Single bit op- 
erations in this system are accomplished by directing dif- 
ferent laser beams to each of the ions, and a fundamental 
gate operation is implemented by exciting the collective 
quantized motion of the ions with lasers, where the ex- 
change of phonons serves as a data bus to transfer quan- 
tum information between the qubits. Prospects of build- 
ing small ion trap quantum computers in the near future 
are good, but the question remains whether stringent re- 
quirements in implementing two bit quantum gates can 
be relaxed: while decoherence of the qubit stored in the 
atomic ground state is not an issue, cooling of ions to the 
vibrational ground state to prepare a pure initial state for 
the collective phonon mode remains a challenge H . The 
question arises whether quantum gates can be performed 
starting from thermal or mixed states of phonon modes. 
The task of implementing quantum computing in "hot" 
systems seems particularly timely in view of the recent 



interest in NMR quantum computing B , where quantum 
states are stored as pseudo-pure states in an ensemble of 
"hot" spins. 

In this Letter we will discuss the implementation of a 
universal two-bit gate between two ions in a linear trap 
with the phonon modes initially in a thermal (mixed) 
state. The novel concept behind the gate operation is 
to implement conditional dynamics based on atom in- 
terferometry of two entangled atoms y|7[]. The two-bit 
gate operation proceeds as follows: ion 2 is kicked left 
or right depending on its internal state with laser light 
(see Fig. 1(a) at i = 0). Thus the ion 1 will experience a 
kick via the Coulomb repulsion conditional to the inter- 
nal state of the first ion. The corresponding wave packet 
will evolve into a superposition of two spatial wavepack- 
ets which are entangled to the internal state of the control 
ion. Provided the spatial splitting of the wave packet of 
the second ion is sufficiently large (at time to in Fig. 1(a)), 
we can manipulate the internal state of the target atom 
depending on its spatial position, i.e depending on the 
state of the first ion, and thus implement a gate opera- 
tion on the qubits. With time these atomic wave packets 
will oscillate in the trap, and with a proper sequence of 
laser pulses this momentum transfered to the two atoms 
can be undone to restore the original motional state of 
the ion (Fig. 1(a) at time tg). The motional state of the 
atom will thus factorize from the internal atomic state 
before and after the gate, independent of being a mixed 
or a pure state. 

In order to analyze the scheme, we consider two ions 
confined in a linear trap [g|. We will assume that the 
motion of the ions is frozen with respect to the y and z 
axis. The Hamiltonian describing this situation is 

H^A + A + vix^) + Vix,) + —^ r, (2) 
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where V denotes the external confining potential along 
the X axis, which we assume to be symmetric V{x) = 
V{—x). We denote by x^ the separation of the ions at 
the equilibrium position. Using the center-off-mass and 
the relative coordinates, 
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and expanding around the equilibrium position Xc — 
cCr = up to second order, we obtain 



H = iJho + Vcoi{Xc, Xr), 
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and Vcor denotes the third and higher-order corrections. 
Here iric — "im and rrir = to/2 are the total mass and the 
reduced mass of the particles, and Vcr the corresponding 
frequencies in the harmonic approximation 
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In the following we will assume that the trapping po- 
tential has been tuned so that v^ = 2i/c. On the other 
hand, in order to use the harmonic approximation to the 
original Hamiltonian we will neglect the corrections Vcor 
in (H). This requires VcoT^g ^ I7 where Vcor denotes the 
typical values of these corrections and tg = iirvc the du- 
ration of the gate. A more quantitative analysis of this 
approximation will be given below. Under these circum- 
stances, the time evolution of the ions motion will be 
strictly periodic, with a period tg] that is. 
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We will denote by |0) and |1) the two internal states 
of atoms which store the quantum information, and by 
1^(0)) = l^int) ® |*mot) the initial state of the ions. Al- 
though we have taken both pure internal |^int) and mo- 
tional states l^mot), our analysis is also valid for mixed 
initial states. The gate is performed in three steps: 

(i) State dependent kick on atom 2: A short laser pulse 
is applied to atom 2 with a k-vector pointing along the 
X direction. The laser intensity is chosen such that the 
internal state of the atom undergoes a flip. The motional 
state of the ions is also changed by this laser beam. We 
write the Hamiltonian for the atom-laser interaction as 
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where 51 is the Rabi fre- 



quency, (T+ = |1)(0| = (<7 )^ and the subscript denotes 
from now on which atom is addressed. The laser is ap- 
plied for a time tias — Tr/ft ^ tg (the so-called strong- 
excitation regime, see Q for details), so that the state 
after this interaction is 



l*'(0)> = (a+e 



ikx2 



-f CTo e 



— ^~ikx2 



)l*(o)>. 



(8) 



According to this expression, if atom 2 is in the state 
|0) it will be transferred to the state |1) and undergo 
a photon kick which will push it to the right. If it is 
initially in |1) it will be transferred to the state |0) and 
undergo a photon kick towards the left [See Fig. 1(a)]. 
Thus, depending on the internal state of the atom it will 
start moving in different directions. Due to the Coulomb 
repulsion, atom 1 will also undergo a different motion 
depending on the initial state of atom 2. That is, the 
motional state of the first atom will also split into two 
wavepackets, ip^' and ip^, moving in different directions 



(right and left). One can easily calculate the evolution 
of the distance between the center of these wavepackets 
using the harmonic Hamiltonian (0) 



d{t) = 2xQrj 
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where xq ~ l/{2mvc)^''^ {Ti = 1) is the ground state size 
of a single ion with the center-of-mass mode frequency 
and 77 — kxQ the corresponding Lamb-Dicke parameter. 
The maximum distance is 



D = rf(io) == 3Vixov/2, 



(10) 



and is produced at tg = 27r/(3i^c)- 

(a) Conditional flip on atom 1: After a time to, a 
laser beam is applied to atom 1, such that it does not 
suffers a kick [pi . The laser is focused on the position x = 
{xe + -D)/2, so that it only overlaps with the wavepacket 
93^, that is, the one that arises if atom 2 was in state 
|0). Adjusting properly the interaction time ii <C ig, it 
induces a rotation |0)i <-> |l)i. The state after this laser 
interaction will be 
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where a'f = a^ + a^ . After this interaction, the atoms 
continue their evolution with the free harmonic oscillator 
Hamiltonian. 

(Hi) State dependent kick of atom 2: After a time tg—to 
a short laser pulse is applied to atom 2 with a k-vector 
pointing along the x direction as in the step (i) for a time 
tias — Tr/fl. Assuming again tias '^ tg, we can write the 
state after this interaction as 

l*(iff)> = {cr^e''''=^ + (72"e-*'=^2)e-*-f^(*»-*")|*(to)>. (12) 
Using (M) one can easily check that 
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which coincides with the fundamental two-bit quantum 

i). 

Under ideal conditions the above steps allow to per- 
form a two-bit quantum gate. According to (^3|) the 
operators acting on the motional state cancel out, and 
therefore the gate can be carried out independent of the 
motional state, regardless of whether it is pure or mixed. 
In a real situation, there will be restrictions to accomplish 
the gate with high fidelity. The main sources of errors 
will be related to: (a) the finite size of the wavepackets 
and their small separations, and (b) the non-harmonic 
corrections to the free Hamiltonian (^. The finite size 
effects can cause several problems during the step (ii): 
(a.l) ip^'^ may overlap at time tg, which will prevent 
us to address one of them alone with the laser beam; 



(a. 2) even if the wavepackets do not overlap, it will be 
hard (if not impossible) to focus a laser beam to such 
small distances; (a. 3) due to the spatial profile of the 
laser beam different positions within the wavepacket will 
see different laser intensities. On the other hand, if the 
wavepackets separate from each other considerably, since 
this seems to be one way to avoid some of the above 
mentioned problems, the anharmonic terms due to Vcor 
may become important. In the following we will address 
all these questions, find the conditions under which these 
problems can be overcome, and present numerical simu- 
lations showing the performance of the scheme in realistic 
setups. 

Let us consider first the problems related to the fi- 
nite size of the atomic wavepackets and their small sep- 
arations in step (ii). We denote by A the size of the 
wavepackets corresponding to atom 1 at time to and by d 
their separation (^) [Fig. 1 (b)] . In order to overcome the 
problem (a.l) it is required D ^ A. The laser profile is 
characterized by the position dependent Rabi frequency 
fl{x) which takes on a maximum value Qq and has a 
width W. In present experiments, due to the impossi- 
bility of focusing laser beams over small distances, this 
width is expected to be much larger than the separation. 
Thus, 



VF > i? > A. 



(14) 



Consequently, the laser beam will affect both wavepack- 
ets ifi^'^, which causes the problem (a. 2). The solution 
to this problem is to select the laser parameters so that 
ipR,L fggj different Rabi frequencies, fl^'^, fulfilling 
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with N integer, the number of complete Rabi cycles. On 
the other hand, we have to make sure that the whole 
wavepacket sees basically the same Rabi frequency, so 
that no information of the internal state is imprinted in 
the motional state, i.e. we have to overcome problem 
(a. 3). According to (n4) this requires that 
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where x = x-^^' is the position of the center of the wave 
packets of atom 1 when interacting with the laser. In 
order to illustrate the above conditions, let us consider 
that the ions are initially in a thermal state, character- 
ized by He, the mean phonon number in the center-of- 
mass mode, and by fir = n^/(2nc + 1) the one in the 
relative motion. For the laser profile we take a Gaus- 
sian n{x) = f^oexp [-(x - 1)^/{2W^)]. We choose the 
equilibrium point of atom 1 to coincide with the steepest 
point of the laser profile, i.e. I = a;e/2 -I- W. The con- 
dition ( |l4| ) can now be expressed as W ^ 3-\/3xo??/2 ^ 



\/nc + nr/2 + 3/4xo taking into account (|lO|). Accord- 
ing to the first inequality, we can expand the Gaussian 
profile il,{x) around x = Xe/2 up to first order. Imposing 
now the second condition (flSh we obtain 



0(Xe/2)^ 



(27V +1/4)^, 
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and W = (4iV -I- 1/2)D. Having this in mind, the first 
condition (n^ can be restated as 



4iV>l, 



> V4rTc-}-2n^ + 3/(3%/3). (18) 



With this, condition (nq) is automatically fulfilled. In 
summary, the laser parameters have to be chosen follow- 
ing the conditions dT^) and (fsl). 

In order to illustrate to what degree the conditions de- 
rived above have to be fulfilled, we have performed some 
numerical calculations. We have considered as an ex- 
ample the potential V{x) = K\x\^/^, which ensures that 
i/j. = 2i'c- We have constructed the evolution accord- 
ing to the harmonic approximation, and calculated the 
averaged fidelity JF and purity V for different tempera- 
tures, i.e. the mean center-of-mass phonon number Uc, 
and values of r]. As shown in Ref. ||l^, these quantities 
characterize the performance of the gate, and the degree 
of decoherence, respectively. In Fig. 2 (a),(b) one can 
clearly see that to obtain a good fidelity for high tem- 
peratures (fie) one has to increase r/ [see Eq. (^]. We 
emphasize that if in a given experiment rj is not large 
enough, one can simply apply a sequence of n pulses, 
from the left and right alternatively, to increase the effec- 
tive displacement of the wavepackets [Q . For example, 
the ^1/2 -^ -05/2 transition of Ca"*" has rj « 0.45, for a 
trap frequency Vc = 2Tr x 50 kHz. In order to obtain an 
effective value oi rj — 7 one should apply a sequence of 
the order of 15 pulses. 

We consider now the second source of errors, namely 
the effects of the anharmonic term l^or that we have 
neglected so far (|^). In order to single out these ef- 
fects, we will assume that the the action of the laser 
beam in step (ii) is ideal. In that case, the errors caused 
by anharmonicities are independent of the internal dy- 
namics, which simplifies the analysis. The fidelity of 
the gate will be is then simply given by the overlap 
jT ^j, = |(^(ig)|*(0))|^ As it should be, if we set ^cor = 
we will have jF(.or = 1 according to (0). Using time- 
dependent perturbation theory we find J^cor = 1~ (AV^)^, 
where 



V 
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(19) 



and (Ay)2 = (*(0)|y2|^(0)) - {^{0)\V\^>{0)f. This fi- 
delity can be evaluated analytically for an initial thermal 
state. In Fig. 2 (c) we have plotted jFcor as a function 
of T] and He- Comparing with Fig. 2 (a),(b), we see that 



the errors produced by anharmonicities can be neglected 
with respect to the ones due to the finite size effects of the 
wavepackets, at least for the values represented in these 
plots. For larger values of rj, however, the anharmonicity 
must be taken into account. 

We have not considered the effects of decoherence dur- 
ing the gate operations, such as thermalization of the 
phonon modes, spontaneous emission etc.. The funda- 
mental limits, however, will allow one to perform many 
gate operations during the decoherence time [p 12|. 

In summary, we have shown how to implement two-bit 
quantum gates between two ions in a linear trap at non- 
zero temperature. The scheme can be easily generalized 
to the case of three ions. In contrast to the ion trap pro- 
posal in ^ the present scheme is not easily scaled up to 
more than three qubits, since one should tune the trap 
potentials so that the frequencies of the different eigen- 
modes become commensurable. We expect the present 
proposal to be of interest in application of quantum logic 
with two and three qubits: examples are fundamental 
experiments involving particle entanglement, demonstra- 
tion of error correction schemes H , and in quantum com- 
munication |13 IJ]. 
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This can be accomplished using two co-propagating 
beams in a Raman transition, or using a laser beam prop- 
agating in a direction perpendicular to the axis x for for- 
bidden transitions. 
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